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Abstract 

| We show that in a large class of stochastic volatility models with additional skew-functions 

(local-stochastic volatility models) the tails of the cumulative distribution of the log-returns 
behave as exp(— c\y\), where c is a positive constant depending on time and on model parameters. 
We obtain this estimate proving a stronger result: using some estimates for the probability that 
| Ito processes remain around a deterministic curve from [3], we lower bound the probability 

that the couple (X, V) remains around a two-dimensional curve up to a given maturity, X 
being the log-return process and V its instantaneous variance. Then we find the optimal curve 
leading to the bounds on the terminal cdf. The method we rely on does not require inversion of 
characteristic functions but works for general coefficients of the underlying SDE (in particular, 
no affine structure is needed). Even though the involved constants are less sharp than the ones 
derived for stochastic volatility models with a particular structure ( [U HZl [9] ) , our lower bounds 
4-H , entail moment explosion, thus implying that Black-Scholes implied volatility always displays 

^-i wings in the considered class of models. In a second part of this paper, using Malliavin calculus 

techniques, we show that an analogous estimate holds for the density of the log-returns as well. 
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^ 1 Introduction 

We consider the following class of diffusions: 
>< : dX t = -\v(t,X t ) 2 f(Vt) a dt + n{t,X t )f(V t )dW? 

(1.1) 

dV t = /3(t,V t )dt + a(t,V t )dW 2 , 

where W 1 and W 2 are two correlated Brownian motions on some filtered probability space (f2, J 7 , 
(•7"t)t>0: 1 )- The function / usually being positive, the couple (X, /(V)) lives in R x R+: when X 
models the logarithm of the forward price of an asset and V its instantaneous variance, Eq. (|1.1[) 
defines a so-called local-stochastic volatility model (LSV). The function rj is the local volatility (or 
skew) function; the autonomous process V is the stochastic volatility. Local-stochastic volatility 
models embed stochastic volatility models (when 77 = 1) and have been intensively studied by 
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the financial community in these last years, in particular when the appearance of derivatives 
whose value depends on the dynamics of the implied volatility demanded the introduction of 
more elaborate models. Some authors have focused on the problem of how to design an efficient 
calibration strategy of such a model to the market smile, as Lipton [H] or Henry-Labordere [12 
[13] , others have given a particular attention to the small-time asymptotics of implied volatility, 
as Forde & Jacquier in [5]. Here we do not focus directly on the problem of the calibration 
of ([l.ip to market data but we rather take up on the issue of giving asymptotic estimates of 
the cumulative distribution and (when existing) of the density of the log forward price X. This 
problem is (clearly) related to the estimation of option prices and the calibration of the model, 
in a precise sense that is to be cleared hereafter. 



1.1 Related work 

Let us precise the context and motivation of our study. Setting F t = Foe * , Fo > 0, then 
F = (F t :t > 0) satisfies 

F t = [ F s r)(s, X s )f(V s )dWg, 



hence F is a positive ltd local martingale. Then, a simple application of Fatou's Lemma allows to 
show that F is actually an integrable supermartingale. Nevertheless, is is well known that within 
the class of models some complications may arise: F may fail to be a true martingale (cf. 

Sin [23], Jourdain [15], Lions & Musiela [20]) and the moments of Ft of order p > 1 may become 
infinite. The latter phenomenon has been attentively studied by several authors (cf. again |20| 
or Andersen & Piterbarg pQ and Keller- Ressel [T7]): if on the one hand the lack of moment 
stability calls for additional care (if, for example, one wants to manipulate variances), on the 
other hand the explosion of moments plays a crucial role in the investigation of the shape of 
the implied volatility surface. Let us recall the basics results concerning implied volatility and 
moment explosion briefly. The model-implied volatility of a European call option with time to 
maturity T and strike Foe k is the unique non-negative solution a(T, k) to the equation 

E{(F Q e x T - F e fc )+] = C BS (k,T,a(T,k)) (1.2) 

where CBs(k,T,<r) is the price of a Black-Scholes call option of strike Foe k , maturity T and 
volatility a. Market-implied volatility cr m arket(T,k) is obtained in the same way, replacing the 
lhs in (|1.2p with the market option prices for different values of k and T. Moreover, let the critical 
exponents p^(X) and qj~{X) of c Xt be given by 

p* T (X) = sup{p > 1 : E[e pXT ] < oo}, 1t{ X ) = su p{<7 > : E[e" 9Xr ] < oo}. 

Lee's moment formula |19| relies the critical exponents and the asymptotic slopes of implied 
variance (the square of implied volatility) in the following way: 

hmsup Tt7( ^ fc)2 = cp(p* T (X)-l), lira sup Tg( ^' fc)2 = y(gr(*)), (1-3) 

k— ^oo ™ k— > — oo ™ 



where <p(x) = 2 — 4(v x 2 + x ~ x),ip(oo) = 0. The relations in (|1.3[) are proved by Lee in a 
complete model-independent framework. We mention here that Benaim and Friz sharpened 
Lee's result: under some technical conditions, they relate the left hand sides in (|1.3p directly to the 
tail asymptotics of the cumulative distribution function of Xt, giving sufficient conditions for the 
limsup to be a true limit. As pointed out by Lee, Eq. (|1.3j) is useful for model selection purposes: 
since the market-implied variance smiles usually display "wings" (i.e. k — > cr mar k e t(T, k) 2 has left 
and right asymptotes), so has to do the model- implied volatility, hence the exponential moments 
of the underlying process X must explode (otherwise p^(X) = q^(X) = oo and (|1.3() tells that 
the implied volatility flattens for large values of |fc|). Pushing things further, the moment formula 
can help model calibration: the values of the slopes at the left hand sides of (|1.3|) can be observed 
on market data for values of \k\ large enough; if on the other hand the critical exponents are 
known functions of the model parameters, the use of (|1.3[) can provide reasonable initial guesses of 
parameters values. This explains the major interest of some authors in the explicit computation 
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of critical exponents, as in [T] or [T7] for some classes of stochastic volatility models: as a result, 
the critical exponents are not available in closed form but can be straightforwardly obtained 
solving (numerically) a simple equation. 

Of course moment explosion is linked to the tail behaviour of the distribution of Xt- More 
precisely, if the law of Xt admits a density and this density behaves as e _c ' y ' for \y\ — > oo for some 
constant c > 0, then positive and negative exponential moments of X t of order p will explode 
for p > c. Dragulescu & Yakovenko [7] first showed that the density of the log-price does behave 
as e c l ^ I in the Heston stochastic volatility model [H] (which is obtained from (jl.ip setting 
r](t,x) = 1, f(v) = y/v, f3(t,v) = k(0 — v), <r(t,v) = a), exploiting the analytical computations 
that can be carried for the characteristic function of Xt- The work of [7] on the Heston model 
has been extended and sharpened with the addition of higher-order terms to the leading e _c ' y ', 
at first by Gulisashvili and Stein in the case of zero correlation and subsequently by Friz et 
al. In [5], relying on affine principles applicable to the Heston model, the authors find sharp 
aymptotic estimates for the forward price density and use them to obtain some fine corrections 
to Gatheral's SVI parametric model of implied variance. 

1.2 Assumptions and results 

The main aim of this work is to show that the cumulative distribution of the log forward price 
and, when existing, its density, behave as exp(— c\y\) for large y in the following class of LSV 
models: 

dX t = -^ V (t,X t ) 2 V t dt + V (t,X t )^V t dW t x (1.4) 



dV t = p(t, V t )dt + a(t, V t )^V t dW t v , (1.5) 

obtained from (jTTTJ) setting f(v) = y/v. We consider finite time horizon T > 0, correlated Brow- 
nian motions d(W x , W v )t = pdt and deterministic initial conditions Xq = and Vb > 0. Eq. 
()1.5[) for the variance process is given on the domain M + = [0, oo), i.e. a process V satisfying 
(|1.5[) is such that P(V* £ [0, oo),t <E [0, T]) = 1. This class contains the Heston model and the 
"universal volatility model" (without the jump part) considered in [21) . but is much wider since 
it allows for general coefficients (3, a in the SDE of the variance. While on the one hand we con- 
sider reasonable Lipschitz, boundedness and ellipticity conditions on the coefficients ?y and a (but 
we allow the drift /3 to be any measurable function with sub-linear growth), on the other hand 
we emphasize that the square-root factors in (|1.4j) - (jl.5|) makes it impossible to apply classical 
Malliavin calculus techniques, or any other method requiring diffusions with globally Lipschitz 
coefficients, to estimate the density of the solution. We also remark that we do not have to deal 
here with the (possibly intricate) discussion on the existence and/or uniqueness of solutions to 
fll.4[l - (|1.5JI : our results hold for any couple of processes (X, V) = (X t ,V t ;t € [0, T]) satisfying 
()1.4[) - (|1.5[) . The situation where the diffusion coefficient of Eq. (|1.5[) is replaced by a{t,v)v p for 
a p > (thus embedding the class of models considered by Andersen & Piterbarg in [T]) will be 
the subject of future work. 

Our main tool to estimate the terminal (at time T) cumulative distribution of X is an esti- 
mation involving the whole trajectory of the couple (X, V) up to time T. In [3] Bally, Fernandez 
& Meda provide estimates for the probability that an Ito process remains in a tube of given 
radius around a given deterministic curve, asking some conditions of local Lipschiz-continuity, 
local boundedness and local ellipticity on the coefficients of the SDE. As a result, the probabil- 
ity of staying in the tube is lower bounded by an integral functional of the curve itself, of the 
deterministic radius and the coefficients of the SDE. In the framework of the model (|1.4JI - (|1.5JI . 
we are able to cast this functional in a simple form and then to optimize over the possible curves 
and radii, obtaining a lower bound which is in the good asymptotic range. To state our main 
result, let us introduce the following objects : for y £ R, let y and the one-dimensional curves 
x t ,Vt,Rt, t G [0,T], be given by 

, / n sinh(f/2) ,„ „. 

y= y +V ; flt = . )' > ; 1.6 

smh(i /2) 
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« t = K)(^0W-e- T / 2 0W+e- t/2 ) 2 ; x t = sign(y)(v t -V ); R t = \ V(V Q A l)v t 
where sign(a;) = 1 if x > and sign(a;) = — 1 if x < 0. Our main result is the following estimate: 
¥(\(X t , V t ) - (Kt,St)| < Rt,0 < t < T) > exp(-c T V>(P±) x \y\j (1.7) 

for | y\ large enough, where ip is an explicit function (cf. (|2.4[) in section 2 for the precise expression) 
and ct is a strictly positive constant depending on the model parameters and explicitly on T 
but not on y nor on the correlation parameter p. The curves x . , v. , R. in (|1.6j) are the product of 
the optimization procedure we set up and they appear as the solution to some Eulcr-Lagrangc 
equations (given in section [3~Tjl . We remark that the curve x. ends up at xt = y while the terminal 
radius Rt is proportional to y/\y\: hence, dropping constants and simply writing P(|Xr — y\ < 
^ 1P(I(^t, Vt) — (xt,vt)\ < Rt) and using (|2.6p . we obtain the desired lower bound for the 
terminal distribution (this argument is made rigorous in Corollary [T] in section |2~2"1) . Notice that 
the fact that the "tube" estimate ()1.7[) is given for the couple (X, V) is crucial in our framework: 
indeed, in order to estimate the behaviour of Xt wc need to have a control on the variance Vt 
for all t G [0,T]. 

In a second part of this work, we go some further and, under some additional hypotheses 
on the coefficients of ()1.4j) - (|1.5j) . we discuss the existence of a (possibly non continuous) density 
for the law of Xt and show that the exponential estimate holds for this density as well. This 
last step requires to work out some "small balls" estimates (cf. Proposition 12.11 in section I2.3j) 
and to employ the integration by parts formula of Malliavin calculus. More precisely, we rely on 
a decomposition of Xt as a Gaussian term plus a perturbation, following the idea of Bally & 
Caramellino in [3J: we obtain the desired lower bound performing a somehow delicate operation 
of balance between the leading Gaussian term and the perturbation, involving regularization of 
the non-Lipschitz coefficients in f|l .4|> - ()1 .5[) and some precise estimates of Sobolev norms of a 
diffusion from [S] . Our final estimates on the density px T of Xt reads 

Px T (y) > — exp(-e T ^(p±)\y\) 

1V±t 

for |y| > My, where Mt and are constants depending on model parameters and explicitly on 
T. If the density is not continuous, the inequality holds almost surely (see Theorem 12.21 for the 
precise statement). 

The paper is organized as follows. In section[5]we give our working hypotheses and a detailed 
presentation of the main results. In particular, in !2.1l we prove estimate ()1.7|) and in !2.2l we state 
the corollary for the terminal cdf, the moment explosion and the implied volatility slopes. In 12.31 
we give our results on the density of Xt- Sections [3] and @] are devoted to the proofs of the results 
stated in section [2J Malliavin calculus only appears in section |H while all the tools employed 
in the previous are borrowed from stochastic calculus for Ito processes. Finally, the Appendix 
contains the proofs of the most technical material and a remainder of the main elements of 
Malliavin differential calculus. 

2 Main results 

In this section we give our working hypotheses and a detailed presentation of the main results. 

Let us consider the class of models (|1.4j) - ()1.5[) . For the ease of computations, wc dccorrclatc 
the driving Brownian motions in the usual way and rewrite (|1.4[) - (|1.5[) as 

dX t = -l V (t,X t ) 2 Vtdt + T)(t,X t )y/Vt(pdW?+p ± dW?), t<T (2.1) 

dV t =p(t,Vt)dt + a(t,Vt)y/V t dW t 1 , t<T (2.2) 

were (Wl, W 2 ; t < T) is a two-dimensional standard Brownian motion. We consider deterministic 
initial conditions Xq = and Vq > 0, finite time horizon T > 0, p E (—1,1) and we denote 
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p± := y/l — p 2 . Eq. (|2.2[) for the variance process is given on the domain R + = [0, oo), i.e. 
a process V satisfying (|2.2[) is such that P(V t G [0,oo),t G [0,T]) = 1. We assume that the 
coefficients ij,f3 and a in (|2.f |) - (|2.2j) satisfy the following conditions, for some K > f : 

(R) (regularity) rj : [0,T]xl-> 1 and <r : [0, T] x [0, oo) — > R are Lipschitz-continuous functions, 
more precisely 

\ v ( s ,x)-ri(t,y)\<K(\x-y\ + \s-t\) 



\<j(s, v) - a(t, u)\ < K(\v -u\ + \s - t\) 

hold for every (s,t,x,y) G [0, T] x [0, T] x R x R, respectively every (s,t,v,u) G [0,T] x 
[0,T] x [0,oo) x [0,oo). 

(G) (growth) The mesurable function /? : [0, T] x [0, oo) — >• R has sub-linear growth in v, more 
precisely 

|0(t,t>)| <*r(l + t;) 

for every (t,v) G [0,T] x [0, oo). Moreover, there exist constants < r\ < 1 < rj and 
< a < 1 < (7 such that 

?7 < ^(^) x) < rj, g_< ait, v) < a 

hold for every (t,x) G [0,T] x R, respectively G [0,T] x [0,oo). 

Remark 2.1. Despite of the boundedness condition on 77, ct given in (G), obviously none of the 
drift and diffusion coefficients in the system (|2.ip - (|2.2[) is bounded, because of the factors Vt and 

Vv t . 

Remark 2.2. In hypothesis (R), we could replace Lipschitz-continuity with respect to the couple 
(t,x) (resp. (t,v)) with Lipschitz-continuity with respect to the state variable x (resp. v) and 
Holder-continuity of exponent 1 /2 with respect to time, and all the results of sections 12.11 and 
would still hold. 



Remark 2.3. We are not interested here in discussing the existence and/or uniqueness of solu- 
tions to (|2.ip ~ (|2.2p under conditions (R) and (G). All the results we give in this subsection and 
in the following (but not in subsection 12.31 where a new set of hypotheses is considered) indeed 
hold for any couple of processes (X, V) = (X t ,V t ;t G [0, T]) which satisfy ([2~Tj) - (j2~2"|) . 

Notation Sets and filtrations. All along the paper, | ■ | will denote the absolute value for real 
numbers as well as the Euclidean norm for vectors, i.e. \x 

I = x i if x e Rn - We denote B R {x) 

the open ball in R™ of center x and radius R, Br(x) = {y G R n : \y — x\ < R}. Moreover, we 
denote {T\,t > 0) the completion of the filtration generated by W z , i = 1,2, and Ft = J-} V J-f . 
A n denotes the Lebesgue measure on R™. 

Classes of functions, derivatives, norms. C ([0,T]) denotes the class of real functions of [0,T] 
which have uniformly continuous derivative on (0,T). We will make use of the class L(/i,h) 
defined in 0], section 2 : given a fixed time horizon T, p > 1 and h > 0, L(fj,, h) is the class of 
functions / : [0, T] — >• R+ = [0, 00) such that for every t, s G [0,T] with \t — s\ < h one has 

/(*) < 

We denote C 0,fe ([0, T] x R") (resp. C°' k ([0,T) xR)) the class of continuous functions of [0,T] xR" 
which have continuous (resp. bounded continuous) partial derivatives with respect to the second 
variable up to order k. Let Qk be the set of multi-indexes of length k with components in 

{1, ...,n}, 9 fc = {1,.. . ,n} k . For a G Ok and g G C°' fc ([0,T] x R") we denote d a g = B ^ . 

We define the norms 

k 

l3|fe = lvVV sup sup \d a g(t,x)\. 

Constants. For a vector of parameters A = (Ai, . . . , \ v ), we shall denote Ca (resp. Ca(£)) a 
positive constant (resp. function of time) depending on the Aj's but not on any of the other 
existing variables. All constants of such a type may vary from line by line, but always depend 
only on A. 
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2.1 Estimates around a Deterministic Curve 

We fix T > 0. We consider three one-dimensional curves x,v,R of class C 1 ([0, T]) such that: 
Rt > for any t £ [0,T] and x and i> have the same initial values as X and V in (|2.ip - (|2.2j) 
(in particular, xq — 0). We look for a lower bound on the probability that a process (X t , Vt) = 
(X t , Vt, t < T) satisfying (|2.1j) - (|2.2[) stays in the tube of radius Rt around the deterministic curve 
(xt,Vt) up to time T, that is a lower bound on the quantity 

P(\(Xt,Vt) - (xt,v t )\ < R t ,0 <t <T). (2.3) 

To lower bound ()2.3[) we employ the estimate provided in j4], Theorem 1. The main result of this 
section, Theorem 12. 1[ makes use of the triplet of curves Xt,i>t,Rt defined in (|1.6|) . As addressed 
in the Introduction, the choice of these particular curves relies on an optimization problem: they 
indeed appear as the solutions of some Eulcr-Lagrange equations (see section 13. ip . Recall the tp 
from the Introduction: 

^•) = ^( ln (J)+ 1 )' r>°- (2-4) 

Theorem 2.1. Assume conditions (R) and (G) and let (X t ,Vt',0 < t < T) be two processes 
satisfying f2.1]) - (K]fy . Then for every y € R with \y\ large enough, precisely 

\y\ > ^ (l + 2sinh(T/2)) 2 , (2.5) 
and with the curves x.,v., R. defined in U.6}) one has 

P(\(X u Vt) - (x t ,v t )\ < Rt,te [0,T]) > exp(-crV(p±) X \y\). (2.6) 
The constant ct is given by 

c T = c*(l + l) e ^ 2 , (2.7) 

where c* is strictly positive constant depending on the model parameters Vq, K, r], g_, rj, a given in 
(R) and (G) but not on y nor on the correlation parameter p. 

Remark 2.4. Let us discuss the impact of the factor ip(p±) and of the maturity T in the lower 
bound (|2.6p a bit further. It is known that the correlation effects moment explosion in stochastic 
volatility models, a negative correlation bringing - as intuitively clear - a dampening effect (cf. p], 
sections 3 and 4). In a Hcston model, obtained when the variance process in (|2.2| has constant 
parameters and mean-reverting drift, the upper critical moment of e XT tends to infinity when 
p — > — 1 and Xt even becomes a bounded random variable when p = — 1, and the behaviour is 
the opposite when p > 1. The factor ip{p±) = ^(yl — p 2 ) has the expected explosive behaviour 
when p — > — 1, but it symmetrically decrements the rhs of (|2.6j) for p > 0, making the lower 
bound significant in particular for p G ( — 1,0). The small time asymptotics ^ of the constant Ct 
is what expected for a diffusion; on the other hand, the large time dependence e c T makes the 
bound ()2.6j) not directly applicable to study the large-time asymptotics. 

2.2 Lower bounds for Cumulative Distribution Function and Moments 

Theorem 12.11 leads, in particular, to lower bounds on the tails of the complementary cumulative 
distribution function (complementary cdf in short) of Xt, i.e. P(|Xt| > •)• Indeed, on the one 
hand we can simply lower bound the probability to be in the tube at the final "time-slice", 
P(\(Xt,Vt) ~ (xt,vt)\ < Rt) with the probability to stay in the tube up to time T. On the 
other hand, the final time radius Rt in (|1.6[) is - roughly speaking - proportional to \/\y\- Hence, 
when y — > oo (resp. y — > — oo) the infimum (resp. the supremum) of the interval [y — Rt, y + Rt] 
becomes large (resp. small) and this allows to obtain tail estimates for P(|Xt| > •) that are in the 
same asymptotic range as (|2.6[) . This observation is made rigorous in the proof of the following 
Corollary, which is indeed a direct consequence of Theorem 12.11 

Corollary 1. Under the assumptions of Theorem \2.1\ for any y > satisfying \2. 5]) and 

y > 2(y Vl) 2 (l + F ), (2.8) 
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one has 

P(A T > y) A ¥(X T < -y) > cxp(-c T V(p±) x y) (2.9) 
where ct is the constant given in \2. 7[ ). In particular, the critical exponents are finite: 

p* T (X)yq* T (X)<c T ^(px), (2.10) 
hence the implied volatility diplays left and right wings, i.e. 



Ta(T k) 2 

limsup -j > ip{cTip(p±) — 1) > 0, 

k— yoo k 

Ta(T k) 2 
limsup > tp(cTip(p±)) > 0. 

k— >— oo ™ 



(2.11) 



Remark 2.5. As addressed in the Introduction, e Xr is intcgrable for every T > 0. A simple 
application of Markov's inequality shows that, for every y > 0, V(Xt > y) = P(e A " T > e y ) < 
e~ y ¥\e XT \. This is not in contradiction with (|2.9j) . because on the one hand ip(p±) i s greater 
than or equal to one for any value of p± £ (0, 1] (cf. (|2.4p ). and on the other the constant c* in 
(|2.7p is greater than 1, hence ct > 1 for every T > 0, too. 

Remark 2.6. In this paper we are mainly interested in the law of Xt- Estimate (|2.6[) can of 
course be applied to derive the analogous lower bound for the joint law of Xt and Vt- See 
Proposition 12 . 1 1 in the next section for a refined statement in this direction. 



2.3 Lower bounds for the density 

We consider now some stronger regularity conditions on the coefficients of (|2.ip ~ (|2.2p : 

(R') {regularity') (R) and (G) hold and r] G C°' 2 ([0,T] x M), a e C*"- 2 ([0,T] x [0,oo)), /3 € 
C°' 2 ([0,T] x [0,oo))niip([0,T] x [0,oo)) with |r?| 2 V |<t| 2 < K. 

Remark 2.7. Under condition (R'), the system (|2.ip ~ (|2.2p admits a unique strong solution. 
Indeed, the existence of a weak solution (X, V) that satisfies (|4.20j) follows from the continuity 
and subdincarity of the coefficients. Then, pathwisc uniqueness holds for (|2.2p after a theorem 
of uniqueness of Yamada and Watanabe (cf. [TB], Prop. 5.2.13) and weak existence and pathwisc 
uniqueness together imply strong existence ([TB], Cor. 5.3.23). Given the unique solution to (|2.2I) . 
standard arguments allow to prove pathwise uniqueness for (|2.ip . 

We will give a lower bound for the density of the law of Xt under hypothesis (R')- Notice 
first that the law of Xt is absolutely continuous with respect to the Lebesgue measure Ai 
on M. This fact may be proven in (at least) two ways. First, we may look to the law of Xt 
conditional to (W^,t < T). Then Xt appears as a functional of the independent Brownian motion 
(W 2 ,t < T) and, using the Boulcau-Hirsch criterium (cf. [22]), we obtain a density px T (W 1 , x) 
for the conditional law. Then, the law of Xt has the density K[px T (W 1 , x)] = px T (x). A second 
way would be to use the results in [BJ (Theorem 2.2) telling that the couple (Xt,Vt) admits a 
density pt(x, v) on M x (0, oo) (meaning that the law of (Xt, Vt) restricted to K x (0, oo) has the 
density pt{x, v)). This immediately yields the existence of a density px T (x) for the marginal law 
of Xt- Nevertheless, we remark that none of the above approaches guarantee that the density of 
Xt is continuous. 

Before giving an estimate of the density of Xt itself, we need to work out some estimates for 
the probability that Xt stays in a ball of "small" radius. 

Proposition 2.1 (Lower bounds for balls of small radius). Let R^\y) be given by 

(so that R {0) (y) = J\y\, R [1) {y) = 1. R (2) {y) = -4==, ■■■)■ Assume (R') and let (X t ,V t ;0 < 

y/\v I 

t < T) be the unique strong solution to Q2.1\) - li2.%\) . Then, for any y satisfying \2. 5\) and \y\ > 
16 V 2(V Vl) 2 (l + Vo) 7 

F(\(X T ,V T )-(y,\y\+Vo)\ <R (j) (y)) >exp(-U + l)dT^(p±) x |»|). (2.12) 
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The constant oIt is given by dx = 2c* + l^e' c +1 ) T , c* being the constant in \2. 7| ). 

Remark 2.8. By taking j large enough, the radius R^\y) can be made arbitrarily small. Then 
we would like to make use of (|2.12p to obtain a lower bound for the density of Xt computed 
at y, but we cannot pass to the limit with j in (|2.12[) because the rhs tends to zero as j — > oo. 
Nevertheless, we can obtain a lower bound for the density using (|2.12j) for finite j and the 
integration by parts formula of Malliavin Calculus. This is what we actually do in order to prove 
the next theorem. 

Here is the main result for this section. 

Theorem 2.2. Assume (R') and let (X t ,Vt;0 < t < T) be the unique strong solution to V2.1\) - 
12.2)) . Then, there exists a strictly positive constant Mt depending on T and on the model pa- 
rameters such that for \\-a.e. y with \y\ > Mt, 

Px T (y) > T exp(-e T ip{p±)\y\) (2.13) 
where er = 136c* + l)e( c * +1 ) T . The inequality V2~W\) is understood in the sense 
[ f{y)px T (y)dy> ^- [ f(y)exp(-e T ip(p±)\y\)dy 

J\y\>M T Mt J\v\>Mt 



for every f G 

Remark 2.9. If the density px T (y) is continuous, then (|2.13p holds for every y with \y\ > Mt- 

We recall that for a Heston model with constant coefficients, the density px T (y) is asymptotic 
to cxp(-c|y|), cf. [313]. 

3 Proof of results in [2J] and [23] 

We start by giving a preliminary result that will be used in the proof of Theorcm l2.ll We consider 
x.,v. in ^([O,^) and R,c,X,j,L : [0, T] -> R+ satisfying 



x = 0; v Q = V ; v' t > 0; 
x' , v 1 , R, c, A, 7, L € L(p, h) 

and we define the stopping time 

r R = t r (X,V) = inf{t < T : \(X u V t ) - (x t ,v t )\ > R t }. 
Moreover, we denote 

b(t, x, v) 

PV(t,x)y/v \ , , _ / p±rj(t,x)^ 



(3.1) 



-^rj{t, x) 2 v 



\ ( prj(t,x)^/v \ , . / p±r](t,3 

a,(t,x,v)=^ j; a 2 (t,x,v)=^ q 



and consider the conditions: 



\b(t,X tATR ,V tATR )\ + Wj(t,X tATR ,V tATR )\ < c t ; (3.2) 

X t I 2 < o-a*(t,X tATR ,V tATR ) < 7i/ 2 ; (3.3) 
[J2 \<r j (.s,X a ,V a )-<r j {t,X t ,V t )\ 2 l {rj 2. a} ] <L 2 t (s-t) (3.4) 



E 

1=1,2 



which correspond to hypothesis (H) in [4]. Then, according to Theorem 1 in [4], the estimate 

„T 

F(\(X t ,Vt)-(x t ,v t )\<R t ,0<t<T)>exp(-Q(ti)(l + J F XtVtR (t)dfy (3.5) 
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holds with the rate function 



.i + ^+«£ +J((4 + ia (i + ^ ) . 



and the constant Q(p) given by 



where 



Q^H-fMn-^, (3.7) 

0a. 7 = inf -; q„ = 8 12 eV 73 - (3.8) 

(We actually denote 0a, 7 the constant p in [3]). The following proposition is the starting point 
to prove Theorem 12. II 

Proposition 3.1. Assume conditions (R) and (G). Let x t ,v tl Rt satisfy h3.1\) and consider a 
process {X u V t ) = (X t , V t ; < t < T) satisfying 12.1 \) -12.2 \) . Let moreover 

Rt<Rv t , te[0,T] (3.9) 

hold for a fixed R € (0, 1) . Then, setting = (K, rj, a, R, Vo), there exist strictly positive constants 
c = ce; L = Lq\ 7 = 7ej A = Ae.r;,^ such that for every < t < s < T the conditions Ii3.2\) - i3.3\) - 
\4\ ) are fulfilled by the curves 

ct = cv t ; L 2 = L T v t ; 
Jt = iv t ; X t = p\Xv t . 

Lt is given by Lt = Le C2T , where C2 is the constant appearing in Lemma j4-2\ ). The curves 
Ct,Lt,^ftjXt belong respectively to L(p,h),L(y/JI,h),L(p,h), L(p,h). 

Proof. In what follows we shall repeatedly apply the inequality ^/v < 1 + v, v > 0. 
(pT2j) : We notice that for every t, x, v € [0, T] x M x [0, 00), 

\b(t,v,x)\ + ^2 WA^ V ^ X )\ 

3=1,2 

< 2^ 2y + K(l + v ) + (P + P±.)vVv + ~a^fv 

< K+ (p + p ± )r] + a + (-rf + K+ (p + p±)rj+a^v 

< c(l + w) 

where the last holds with c = ^rj 2 + K + 2rj + a. Then, employing the condition ()3.9|) on the 
radius and the fact that v t >Vq for any t € [0, T] by (|3.1|) . 

3=1,2 

<c(l + V fATR )<c(l + (v t +R t )) (3.10) 

< c(l + R)(l + v t ) < 2c(l + fljiX-O^ 

^0 

and the last inequality holds since (1 + v) < 2^-v for any v > V . 

(|3.3p : Let aa* -(t, z, w) = 5Zfe=i 2 "fcC^ ^ u )°fc(^ x i w )> J = 1> 2. The condition 

\ t I 2 < aa*(t,X tATR ,V tATR ) < 1th (3.11) 

for the given At, 74 will follow from the computation of the eigenvalues of aa* . Denoting 77 := 
rj(t,x) and a := o~(t,v) for simplicity of notation, we have 



aa* (t,x, v) 



rj 2 v prjav 
pr\av a 2 v 
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hence the smallest, respectively the largest, eigenvalue satisfy 

At (a;, v) = — (^i] 2 v + a 2 v — \J (jfv + a 2 v) 2 — Affa 2 v 2 p± 



2 >jW 2 gV 



' 77 2 W + <7 2 t> 



(3.12) 



7 t (x, u) = — (^i] 2 v + a 2 v + \J (?] 2 v + fT 2 u) 2 + 4?] 2 a 2 v 2 p± 

^2,2^ 1—2 . — 2\ 

<r;t; + (Tu<(?7 + cr Jt; 
Proceeding as before we have 



a 2( ZZ 2 g 2 ) 1 - R 
rj 2 + a 2 \ + R 



v t = p ± X(r],g_,T],a,R)vt; 



j t (X tATR ,V tATR ) < (t + a 2 )(l + R)v t = 1 (Ti,a 7 R)v t ; 
Then (|3.11l) follows with At, -ft as in the statement of the proposition. 

(|3.4[) : Because of assumption (R), for every s,t € [0,T] x [0,T], every x, y £ 1x1 and every 
v, u 6 [0, oo) x [0, oo) we have 



^(SjX)-/;; - r/(t,y)y/u\ < y/uK(\x - y\ + \s - t\) + r)\ y/v - y/u\ 
< y/uK{\x-y\ + \s-t\ 



\v — u\ 



(3.13) 



and 



\a(s, v)\/v — a(t, u)\/u\ < \/uK(\x — y\ + \s — t\) + a\y/v - y/u\ 



< V^K(\x -y\ + \s - t\) + \v - u|. 

2 mm(y/v, \/u) 



(3.14) 



It follows, for every t < s < T, 



E 



Ms, x.) Vv s - f](t, X t ) yJV t \ 2 l {TR > s} 

< 4K 2 E[V t (\X s - X t \ 2 + (s- t) 2 )l {TR > s} ] 



E[ 



< AK 



'■[(l + R)v t C 2 e c ^ 2 
,VnVl 



2mim%V s ; 



t\ L {tr>s}\ 



is-t) + T { s-t)) + ^ w C 2 e^\s-t) 



v t (s - t) 



1(1 -R)V 2 , 

where C 2 is the constant considered in Lemma 14.21 Analogously, 
\a(s,V s )y/V s - a(t,V t )VV t \ 2 l {TR > s} 



2 



< 4K 2 ( (1 + R)v t C 2 e C2T (s-t) + T{s - t) 



a 2 



2(1 - R)v t 



< C 2 [S(1 + R)K 



2 VbVl 



7f 2 



V (1 - R)V 2 



v t {s - t). 



Estimate (|3.4p then follows from the two previous inequalities and the expression of cti, a 2 - 
The last statement on the curves q, it, 7t, At follows from the fact that the function af p belongs 
to L(/j p , h) if / belongs to h), p > and a is a positive constant. □ 

Basically, what Theorem 12 . 1 1 docs is to compute the right hand side of (|3.5|) on a particular 
curve satisfying conditions (|3.2[) - (|3.3p - (|3.4p . so that (|3.5I) translates into the explicit lower bound 
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(|2.6[) . The choice of the deterministic curve (xt, vt) considered in Theorem 12. II is of course moti- 
vated by the form of the rate function (|3.6[) . More precisely, consider any xt,vt,Rt that satisfy 
(j3.1[) . Then, by Proposition 13. 11 the estimate (|3.5[) holds with 

W(t) = 1 + MM + 2(c ^ + LTVt) (_L_ + _L) , (3 ,5) 

0A7 = inf ^ = d^ inf ^ = d^ ( 3.i6) 



and Q(fi) as given in p. 71 



3.1 A Lagrangian minimization problem 

We start from the simple observation that maximizing the lower bound in ([3.5)1 is equivalent to 
minimizing the exponent Q(/x)(l + J Q F X}V _ft(t)dt). Due to the presence of the competing terms 
— + -g? and {x' t ) 2 + (v' t ) 2 in F x , Vj r, we make the choice 



Rt^^VVQ^t (3.17) 

so that R 2 is proportional to Vt, and consider curves Xt,v t such that \x' t \ = \v' t \, precisely 

x t = sign(y)(v t -V ), ie[0,T]. (3.18) 

Equations (|3. 17|> and (|3.18p define R t and xt given v t , as happens for (jTTBJ) . We remark that the 
radius in ([3.17P satisfies the requirement R t < ivt of Proposition 13.11 Moreover, if the arrival 
point xt = y of the curve xt is given, the same will be for vt- We define the "shifted" arrival 
point y = vt setting 

y = \y\ + v . (3.19) 

After ([337]) - ([535]) . the rate function ([3~TS]) reduces to 

F v (t) = i + ^K) 2 + 2(c 2 ,? + L TVt )(±- x +1)1 (3.20) 

which is a function of the curve «t only. Since we want to upper bound F v , we can get rid of all 
the constants and just keep the explicit dependence with respect to the curve vt- defining 

r T = i v2 /ta c \ Lt „ - , (3.21) 

(y Al)piA' 



have 



h V v t 

and the constant T carries the explicit dependence w.r.t the model parameter p±. The strategy 
we shall follow is to consider 

C(v u v' t ) = ^-+v t (3.22) 
v t 

and to look for the solution of the minimization problem 



mm 



C(v t ,v' t ) (3-23) 



with the constraints 

^0=^0; v T = y. (3.24) 

The problem ([3.23p ~ ([3.24p is the classical minimization problem in Calculus of Variations for 
Lagrangian systems: a stationary point for the integral functional in p. 231 ) is given by the 
solution of the Eulcr-Lagrangc equation 

d dC , dC , 
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under the constraints (|3.24[) . Some simple calculations yield the Eulcr-Lagrange equation asso- 
ciated to the Lagrangian (|3.22[) : this equation reads 

— = 7T- + T7- 3 - 25 

v t 2v t 2v t 

A closer look to Eq. (|3.25[) reveals that it can be converted into a linear second order ODE - 
hence explicitly solved - with the change of variables 



which indeed converts (|3.25[) into 



7«t = 0, (3.27) 



4 



now with the constraints 



_ i 

V 



u = 1; u T = (^J . (3.28) 



The explicit solution to ()3.27|) - ()3.28|) is easily found to be 

y_\ \ sinhft/2) -r/2 sinh(f/2) _ t/2 



\V J sinh(T/2) sinh(T/2) + ' [ ' 



The curve Vt defined in (|1.6|) corresponds to the one given by (|3.26p and p.29[) . What Theorem 
l2.1l does. then, is to pick up this particular curve, to check for which values of u, h and y the curve 
v' t belongs to L(/j,, h) and satisfies v' t > from (|3.ip . hence to estimate the integral functional in 
^Mi- 
Proof of Theorem\MJ\ Step 1. We show that v' t > 0, t £ [0,T] and v' £ L(A,h) with h = 
i 

(v0 2 tanh(T/2), if y satisfies (|2.5p . Taking advantage of the notation introduced in Q3.26p . we 

v' t = 2Vbwtw'f 



have 

and a simple calculation yields 



y_y _ p -T/2\ coshft/2) _ 1 _ t/2 



«,/ e "'"JscE^j-?""- (3 - 30) 

We remark that w t > for every t £ [0, T] as soon as y > X , hence by (|3.27p u" > 0, too, and 
consequently u' t is an increasing function. Using the expression of u' given by (|3.30[) , it is easy 
to verify that (|2.5p implies u' t > u' Q > j > 0. Now, we simply observe that u £ L(2, Hu'H^ 1 ): 
indeed, for every s,t £ [0,T] such that |s — 1| < 



-l 

I CX) ' 



U s < U t + \\u \\oo\s — t\ < U t + 1 < 2Ut 

where the last holds because u t > uq = 1. Analogously, v! £ L(2, \ \u\\^) because 
u' s <i4 + ||«"||oo|*-t| <t4 + (l-o) 2 ||u|| 00 |«-t| <2u' t 

holds if \s — 1\ < llull^, 1 , employing in the last step the fact that u' t > u' > \. Because u t and u' t 
are increasing, we have 



— , 1 

y \ 2 i 



& \ 1 1 / 1 1 ^ 

Ut = I "TT ) i II" Hoc = u T < 



,y / 2tanh(T/2) 

Observing that tanh(T/2) < 1, we conclude that both u and u' belong to the class L(2, h) with 

i 

h = (lb) 2 tanh(T/2). The fact that v' t £ L(4,h) for the same h now follows from the property 



cfg £ L([if x Li g ,hf A h g ) if / £ L(fj,f,hf), g £ L(/z ff ,/i ff ) and c is a constant. 
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Step2. We estimate the integral functional at the right hand side of (|3.5[) . 
By Proposition 13.11 and the computations at the beginning of the current section, we know that 
the rate function F~~^ is upper bounded by Fj, defined in (|3.20[) . more precisely F- ~ ^ < 

\ + — V Vtj ■ Making once again use of u defined in (|3.26l) . we have {v' t ) 2 = AVQU 2 (u' t ) 2 , 



hence 



r^(^<| J (i + r T (M + , t ))^ 



<| + r T / T f 4 y ^Ml! + 

T 

77 



< Z + 4T T V I ((u' t ) 2 + u 2 t )dt 



< 



(i) 4 ss^72) +4r ^"r«"» i+ «->' a 



T 

tanh(T/2) 

and we just have to integrate the expressions for u t and itj over [0, T]. Since we are interested in 
an upper bound for the integral, we simplify the computations using 



( l\5(sinlf2) \ , 1 (V \h cosh(t/2) 

: \ y / V sinh(T/2) + J ' 4 " 2 \V ) sinh(T/2) ' 



Hence, setting 



V S inh(t/2) n» (2) _1 1 [ T ^umM. 



c - = J + X J ^ - iMfw 7o cosh(t/2) dt 



T 



■ 4y (4 1) + 4 2) )) 



,tanh(T/2) 

we obtain that, if ?/ satisfies (|2.5|1 so that in particular (^r) 1 ' 2 < ^ < 2|j/| 



/ F £fi ^)dt<c T r T N = c T r T |y| 

Jo 



We remark that have ci^ < fJ 4dt < AT, cl 2) < 1 - ? 2 < ^ + T and f JL,^ < T + 1, hence 
ct < 4(20Vb + 1)(A + 7 1 ) for a positive constant c depending on Vq. On the other hand, recalling 

the expression of IV from (13.211) . we have Tt < -^-e C2T for a positive constant T > 1 depending 

p± 

on V ,k,T], a, r},W but not on p± or T, hence c T T T < 4(20V" + l)^-(^+T)e C2T2 < -^(^ + l)e c * T2 
with c* > 1. Now, by p,16[) . the constant Q(/j,) in (|3.7|) is given by 

= T4-TI ln - 



with g = 8 12 e 2 4 73 . Eventually multiplying the constant c* by 3^ In we conclude that 
exp 



(-qo*)(i+ y o /•. ,,,,■:' w//)) < 

cxp(-2c*(i + l) e c * T2 -l(ln(-^)+l) x |y|) (3.31) 

for every y satisfying (|2.5j) . Using (jHHJ) and the definition of i\) in (|2.4j) . the proof is completed. 

□ 



We now prove Corollary Q] 
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Proof of Corollary^ We consider y € R* with \y\ > (1 — Vq)/2 and we now define y = 2\y\ + Vq 
and consider x t , Vt,Rt as in Q1.6p . We remark that Rt = \\fVa A ry2|y] + Vq < | if \y\ is larger 
that the larger root of \y\ 2 V~o — 2\y\ — Vq. This holds in particular if \y\ > 2(Vq V 1)(1 + Vq) as in 
JUHJ. If y > 0, we write 

P(*t > if) > P(|X T - 2y\ < |) >P(|X T -2y| < Rr) >F(\(X T ,V T ) - (2y,V)\ < Rt) 

and 

P(^t < -y) > P(|X T - (-2y)| < y) > P(|X T - (-2y)| < JZr) 

>p(|(x T ,VT)-(-2»,y)| <-Rt) 

and in both cases the last term is larger than F(\(X t ,Vt) — (xt,vt)\ < Rt,t £ [Q,T]). Theorem 
Othen yields estimate ([23]) . To prove ([2~TU]) . we shall first show that if E[e pXT ] < oo, then 
V < CTip(p±)- Indeed, it is sufficient to observe that if E[e pA " T ] = C < oo, p > 0, then ¥(Xt > 
y) < Ce~ pxy for all y > by Markov' inequality: 

P(X T > y) = P(e pA " T > e py ) < e - pxy E[e pXT }. (3.32) 

Since (|2.9j) and (|3.32p hold simultaneously for all y from a certain range on, clearly this implies 
P < CTip(p±)- With the same argument and using the estimate for F(Xt < —y), one shows that 
if E[e~ 9A " T ] < oo, q > 0, then q < CTip(p±). Finally, the estimate (|2.1ip on the implied volatility 
is a direct consequence of moment formula (|! -3p and of (|2.10p . recalling that the function ip is 
decreasing. □ 



4 Proof of results in 2.3 



We introduce some compact notation that will be used throughout this section. For t, s with 
< t < s < T and G R x [0,oo) we denote , V^ 1 ; i < u < s) the solution of 

(|2.ip - (|2.2p on [t,s] with initial conditions X t = X\ and Vt = v\. We denote Y*' Xl,Vl the couple 
(X^,^) and 



^' X2 = Xl + X2 Xl (u-t), ue[t,s] 

s — t 

C'" 2 =vi + ^jp t -(u-t), «eM 



(4.1) 



the line segments between (i, xi),(s, #2) and (t, t>i),(s, W2) respectively. For y 7^ and a couple of 
radii R2 with < i?2 < -Ri < |y|, we define 



^,(2/,^,^)= inf P^ 1 / 1 ^,^)). 

(aJi.uiJefiRiCv.lyH-Vo) 

Moreover, we set 

O \T] 0~ 

eo = -7=^ A 1 (4.2) 

with eo = 1 if p = 0, and 

^iJ^f; '= p fe ws ili) (43) 

where (b u ,u > 0) is a standard Brownian motion under P. The following lemma provides some 
estimates that will be used in the proof of Proposition 12.11 and Theorem [ 



Lemma 4.1. Let y € K with \y\ > 16 and Ri,R 2 with < R 2 < Ri < yjyj- Assume (R'). 
Then, for any < i < s < T , 

pUy, Ri,R 2 ) > exp(-c T V;(p±) (j^^ + - *)) ) (4-4) 
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where ct is the constant defined in {2. 7| ). Moreover, ify > 0, for any t > and any < 5 < ^ 



At 



we have 

inf 

" e - B «ov^/2(y) 



°{\V*' v -y\ <e Q ^6y~,t-6 <s<t; 



t-s 



{a{uXn-^{t-^ s ))\f^dWl\<e ^)>\q. (4.5) 



The proof of this lemma is not particularly enlightening for the rest of our study, hence we 
postpone it to Appendix 14.21 Here we give the proof of Proposition 12.11 

Proof (of Provosition\2.1\). Step 1. We consider Ri,R 2 with < Ri < R 2 < \J\y\ and ^ < t < 

s < T. We have {Y s € Br 2 (y, \y\+V )} D {Y t G B Rl (y,\y\+V Q )}n{Y u G B R2 (x*^, vl tM+Vo ), t < 
u < s}. Hence, applying Markov property for the process Y 



1 {Y t eB Rl (y,\y\+V )} E 

1 {neB Rl ( !/ ,|2/|+vb)} IE 



\Y S G B R2 (y, \y\ + V )) > P({Y t G B Rl (y, \y\ + V )} fl {Y u G B Ra (x*" u , v?>M +v °), t < u < a}) 

1 {Y u eB R2 (xXf\vZ'M+ v o ) ,t<u< s }\ :Ft _ 

>P(rteflfl 1 (y,|y| + y ))xpt,.(»,i2i,i22)- 
Step 2. We define the time step 

T 



(4.6) 



Applying Lemma ()4.1[) . for any j > 1 we have 



i>i. 



;< inf ^(y,^-^),^^)) > ex P (-c^W)(^M! + 



exp 



(-c t V>(p_l) 



V (i?«(2/)) 2 JV 

y 2 - J y 2 T 



T\y\i-i \y\^2\y\o 
exp(-2c T ^(p±)(i+T)|y| 



(4.7) 

On the other hand, R = \y/{V A l)(|y| + V Q ) < ±y/\y\ + V Q < < R {0) (y)- Applying 

Theorem (|2.ip on the interval [0,f], we have 

inf P(F t G flawed/, |]/| + Vo)) > inf P(Y t G B R (y, \y\ + V )) 

- T^< T cxp ( _c *(i + l )^p±)\y\) (4-8) 

>exp(-2c T (i + l)v(p±)M)- 

Step 3. We fix j € N* and define 

k 

t> k = T-Y, S i-h+i, 0<k<j 7 



h=l 



so that to — T and t J k _ 1 — t J k = 8j-k+i for 1 < fc < j. Moreover, since Ylk=i — £fe=i "T — 
T^pT < f , we have t J k > ^ for all <Kj. Repeatedly applying (|4~E| and (|4~7) . we 
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get 

P(Y T G B RU) (y, \y\ + V )) = P(Y to 6 B rU) (y, \y\ + V )) 

> p(Y tj g s^o)^, \y\ + v )) f[p tltl Jy,R u - k) (y),R u - k+1] (y)) 

fe=i 

> P(F t . G Bh(o)(i/, ||/| + V )) x cxp(-2jc T (i + r)v(p±)|»|) 

> exp(-2(j + l)c T (i + T V 

and in the last step we have applied (|4.8|l . Using the expression for the the constant ct given in 
1(27)1 . we have c T (y+TVl) < c*(i + l)e c * T2 (i + rvi) < 2c* + l) e ( c *+ 1 ) T2 and (l2~T21) 
is proved. □ 



Let us go back Theorem 12.21 To lower bound the density of Xt we follow the approach of [3J , 
section 5. The idea is to treat Xt as a random variable of the form 

F = x + G + R, 

where x G K, i? G D 2 -°° and G is a Wiener integral G = E,=i,2 J T h A t ) dw t \ witn : 
[0, oo ) — > R deterministic. Here D 2 '°° denotes the space of the random variables which are two 
times Malliavin differentiable in LP for every p > 2 (we refer to [22] for the notation and for 
a general presentation of Malliavin calculus; see also Appendix 14.31 for a reminder of the main 
elements of this theory) . Remark that G is a centered Gaussian random variable with variance 
A = J^j Jo " hj(t) 2 dt > 0. Let g^{y) = ^t-A cx p(~2/\) denote the density of G and ||-R||2,p the 
stochastic Sobolev norm of R of order two. Our starting point is the following result due to Bally 
and Caramcllino in [3], which we restate here in a form suitable for our purposes. 

Proposition 4.1 (Proposition 8 in (3j). If the law of F has a density pf , then for any f £ G&(R) 
one has 

" f(y)p F {y) > I f(y){9A(y -x)- e(A, R))dy, (4.9) 
r Jr. 

with 

e(A,R) = ^(l + \\R A \\ 2 ^f\\R A \\ 2 ^ 



where R A = R/y/A and C*,q*,l* are universal constants. 

Proof. Using point i) of Proposition 8 in [3J , we know that there exists a probability measure P 
on (fl, J 7 ) such that jjjp < 1 and the law of F under P is absolutely continuous with respect to 
the Lebesgue measure. Again according to [3J, the associated density p F satisfies 

sup \p F (y) - 9A(y ~x)\< e(A, R) 

2/6R 

for the given e(A,R). (We refer to [3J for the explicit construction of the probability P). Then, 
for any / G G&(R) we have 

f(y) PF {y) = E[f(F)] > e[/(F)^] > [ f(y){g A (y - x) - e(A, R))dy 



R 



R 



which proves ()4.9p . 

Remark 4.1. If the density pf is continuous, then (|4.9[) implies pf(v) > 9a(v — %) — e (A, R) 
for all y G R. 
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Remark 4.2. Wc shall use conditional calculus in order to prove Theorem l2.2l : in particular, we 
will work with Malliavin derivatives only with respect to the Brownian noise Wt , t £ [T — 5, T] , 
and consider conditional expectations with respect to J~t-s, for a 5 < T. This allows us to gain a 
free parameter S in (|4.9p that we can eventually optimize, and this feature turns out to be crucial 
in our analysis (cf. Propositions 14.21 and 14.31 hereafter) . The use of conditional Malliavin calculus 
in order to derive lower bounds for the density of a random variable is not new and has been 
employed by, among others, [T5], [5] and [TU]. In our framework, we face some supplementary 
difficulties. Let us point them out: first, to estimate the marginal density of Xt we have to 
separately estimate the whole path of the stochastic volatility V up to time T. This was the 
motivation of estimate (|4.5j) . Second, in order to manipulate the Sobolev norms of R we need all 
the involved random variables to be smooth in Malliavin sense, but this is not guaranteed in our 
framework due to the presence of the non-Lipschitz square-root coefficients in (|2.1[) and (|2.2I) . 
This is why we introduce a regularization of the coefficients of the SDE, as we do hereafter. 



Let us implement what stated in Remark 14.21 We consider the case of positive y in Theorem 
2.21 the case of negative y is proven in the analogous manner. We assume y > 2 and introduce 
two parameters S > and I £ N such that: 

So 11 
y y L 2 



S < ^; - < ^oP±WyS (4-10) 



and eo , Sq as defined in ()4.2|) and (|4.3p . We remark that for such a value of 5 we have eo \fb~y < 
eoV^o < 1- Then, we consider a truncation function ip £ C£°(R,R) such that ip(x) = x for 
| — j/ 1 < 1, ip(x) = y — | for x < y — 2 and i/j( x ) — y + § f° r x > y + 2. if) can be defined in such 

a way that j^lo < %J + I < 2y and X^=i su PzeR IV^H^)! — 2 ' 2 ' ■ We define the sets 

As t i(x, v) = {\x T _8 -v\<- v \Vr-8 - (y + v )\< e -^VyS} 

W 2 

and 

A S (V) = {\V s -y\< eo^j.T- 6<s<T; 

(o-(«, K) - <t(T - 6, V T -s))VV u dWi\ < eoV^}, 



,7' 
/T-5 

and denote 

A s ,i = A 5 j{X, V) := A Stl (X, V) n ^j(V). 
Finally, we consider (X t , Ft; T — # < i < T) the (unique strong) solution to the equation 

r 



X t =X T _ 5 -\( v (s,X s ) 2 T/j(V s )ds+ [ r](s,X s )M(V s )(pdW^ + p±dWf), 

■ 'T-5 Jt-S 

(4.11) 



V t = Vt-8 



+ f P(s,V s )ds+ f a{ Sl V s )^{V s )dWl 



T-5 JT-5 



We remark that on the set A$j, ip(Vt) = Vt for all t £ [T — 6, T}. Hence, since pathwise uniqueness 
holds for (|iTTj) . we have (Xt, Vt)( w ) = (Xt,V t )(oj) for (t,ujj_£ [T -_S,T] x A S j and in particular 
A s ,i = A s ,i(X,V) C A s ,i(X,V). Under hypothesis (R'), X t and V t belong to the space P 2 ' p 
associated' to (W?, W?),'t £ [T - S, T^for all p > 1. 
We decompose the random variable in the following way: 

Xt = Go + G + R, 
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where 



G = X T -s + pv(T - 6, X T -s 

r-T 



T 
T-S 



ip{V t )dWt 



G = P± r 1 (T-5,X T _ 5 ) 

i-T 



T-S 



iKv t )dw? 



R 



(4.12) 



V (t,X t ) 2 ip(V t )dt 



T-S 



+ / fa(t, X t ) - r](T - 5, X T -s))y/1>{V t )(pdW} + p±dWt). 

J T—S 

Conditional to Tt-s V T\, the random variable G is a centered Gaussian with variance / = 
p 2 ± r](T — 5, X T -s) 2 $T_ s il){Vt)dt. By the definition of ip, we have / > p\t} 2 J^_ 5 (y — \)dt > 



\p 2 ± r i 2 y&- Similarly, we can see that an upper bound for I is given by 2p 2 L rj 2 y8, hence 



with 



A < I = Var(G\F T -s V 7£) < aA 
1 



(4.13) 



a = 4 



7] 2 



Using (|4.13[) and Lemma ()4.1j) . we can prove the following statement (which is the analogous of 
Lemma 5 in [5]): 

Proposition 4.2. Let g(-\Fx-s V F^) denote the density of G conditional to Ft-s V F^. Then, 
for any y, 5, 1 satisfying 

1 



g(y-G \F T _ s V J*) > 



p±?]ey/4ir5y 



on the set As 



(4.14) 



Proof. Recall that I = p±T}(T - S,X T -s) 2 J^_ s ^{V t )dt. Moreover, let us set J = pr/(T - 



(4.15) 



S, X T -s) St_ s yfi>(V t )dWi (so that G = X T _ S + J). Then 

g(y - G |F T _ 5 V 7*) = -^=exp(-l(y - (X T _ S + J)) 2 ) ■ 
Since / > A and \y — Xts\ < on Agj, on this set we have 



\y - Xts\ 



< 



< 1 



VI " VA 

where the last inequality holds because of (|4.10l) . Now, using equation (|4.11[) for V, 



<t(T - (5, V T -s) 



T-S 



^{V t )dWl =V T - V T -s - 



T-S 



P(t,V t )dt 



hence, on the set Agj 

rT 



T 

T-S 



ip{Vt)dW t l < -[\Vt-V T -s\ + 
s • £ v 

T 

T-S 



(a(t, V t ) - a(T - 5, V T -s))^{V t )dW^ 



\/3(t,V t )\dt 



T-S 



(a(t, V t ) - a(T - S, V T - S ))J iP{V t )dW} ) 

< - (2e \/yS + K{l + y + e V^)<5 + eo V^) 

< -{2e Q + 2Ky/y~5 + e Q )^y~5< — y/yS 
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and the two last inequality are obtained using K(l + y + eo^/y5)6 < K(l + y+ 1)5 < 2KyS, then 
2K^/yS < 2Kji~Q < e after (pi~3) . The previous estimate yields \ J\ < ^^■y/yS', hence 

|J| |J| 4V2^ 

-7= S —7= S £o S l 

V/ VA 

and the last inequality holds after (|4.2[) . Finally, for the exponential term in (|4. 1 5[) we have 

cx P (-^(y- (Xr-* - J))") > e -5(W) 2 > e - 2 (4.16) 
on the set A 5>1 . Since / < a A = 2p 2 A rjy8, (|4TT5jl yields (|3TI3|l . □ 

The second result we need in order to prove Theorem 12.21 is an estimation of the reminder R. 
Let i?A := R/y/A as in Proposition 14. H 

Proposition 4.3. Let y, 6,1 satisfy J^. lOty . Then, for every p > 1 there exists a positive constant 
c p such that 

\\Ra\\t-s.s,2, p < —V&y^x e CpTP on the set A S r. (4.17) 
P±V 

The constant c p depends also on K but not on the other model parameters. 

Remark 4.3. Similar estimates (with different powers of i5 and y) could be obtained for 
H^?a||t-i5.i5,/£,p, k > 2, under the corresponding regularity assumptions on the coefficients 77,/? 
and a. 



We can now prove Theorem [ 

Proof (of Theorem \2.2\) . We make an explicit choice of 5 and I: 

p 1 77 2 2~ r_5 / 2 

6 = S @T x?T 31 ; ' = 16; with & T = t_ _^ e- 2c » TP (4.18) 

y/ireriC*Cp 

where C* and c p are the constant appearing in, respectively, (|4.9|) and (|4.17[) . Conditions (|4.10[) 
are satisfied as soon as 

y > Mt with Mt = 2/(eoP±'T]^/So®T)- 

We now apply Proposition 14.11 to the law of Xt conditional to T = J~t-s V J- T - Let p^ T {-\F) 
(resp. gOI-? 7 )) denote the conditional density of Xt (resp. G), by Prop. POl we have 

/ f{y)PxM y ) d v^ I f(y)(g(y-Go\T)-e(A,R))dy 

Jy>M T Jy>M T 

>[ f(y)( _ 1 -e(A,R))dy 
J y >M T \pj_r]e^4TTdy I 



with 



e (A, J R) = -^(l + E[||i? A || 2 , 9 ,|J-]) r E[|| J R A || 2 ^|^] 



-^l(l + E[||i? A || T _ 5A2 . 9 ,|J-i]) r E[||i? A || T _o,|J"T] 
V A 



For the given value of 6, \\Ra\\ts s 2 9 * < 1 on the set A$i, hence by Prop. FOl ef A, R) is bounded 

by 



e(A, R) < 5% x 2<* +1 / 2 x y/t^e ' 7 * 
PIT 



on the set A$j. The value of i5 in ()4. 18[) is chosen in such a way that the right hand side in this 
last estimate is smaller than i x — - , , hence 

I f{y)PxM y ) d v^ I . J (v LcJ y= I , - y 15d y 



ly>M T Jy>M- 



T p^ne^/AirSy J v> m t p±r] e \ / ^6 ®T' 
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on the set Agj. 

Let us now estimate the probability of the set Agj. Since -\ < by condition (|4. 10[) . 

then Ats 3 {(Xts, Vr-s) G B±(y,y + Vb)}. Hence, we apply Prop. Hj] for j = 21 + 1 and 
obtain 

P{A T -s)>P{{Xt-5,Vt-s) € B^(iM/ + Vb)) > cxp(-(2? + 2)d T _5V(p±)2/) 

Since 2/ + 2 = 34 and (T ^ )2 < ^, then (21 + 2)d T s < 136c* (-^ + l) e ( c * +1 ) T := e T and 

F(At-s) > exp(— eTip(p±)y)- Applying Markov property for the process V and (|4.5p in Lemma 
14.11 it is easy to see that 

HAsj) > |P(A T _ 5 ) > | exp(-e T ip{ P ±)y) . (4.19) 

Finally, let us denote px T (-\J~) the density of Xt conditional to J- '. We have 

f f(y)px T (y)dy= [ f(y)E{ PXT (y\T)}dy 
Jr Jr 

> f f(y)E[p XT (y\T)l Asi ]dy 
Jr 

= / f(y)E[p TT (y\T)l A Jdy 
Jr 

> f f(y) _ * n -2/ 15 x nAs,i)dy- 
Jr p_Lr?eV47rd BT 

Using estimate P~H?)) . we obtain ([2~T^| . □ 

Conclusion 

We have shown that the left and right tails of the distribution of the log price X decay no 
faster than exponentials in local stochastic volatility models driven by square root diffusions 
- namely, in the model class (|1.4[) - (|1.5[) - no matter how the (possibly time-dependent) skew 
function 77, the volatility drift (3 and volatility of variance a are chosen, provided they satisfy 
some reasonable boundedness and linear-growth conditions - namely, conditions (R) and (G) in 
section [5J Together with the elementary observation that e x is an intcgrable supermartingalc, 
this yields the "sandwich" estimate e~ Cl ^ v < F(X t > y) < e~ C2 ^ v for large values of y. From 
the point of view of the financial modclisation, our estimate has an impact on moment explosion 
and, by Lee's moment formula, it turns into lower bounds on the asymptotic slopes of the implied 
volatility. 

Our result is not limited to fixed-time marginal laws: we have shown that the exponential 
lower bound actually holds for the probability that the whole trajectory of the couple (X, V) 
remains in a "tube" of given deterministic radius around a given deterministic curve for all the 
times up to a given maturity. This means that our main estimate can also be applied to the two- 
dimensional joint distribution of X and V and to study the law of suprema of the components of 
the solution to (|1.4j) - (|1.5|) . Back to the financial level, this can eventually lead to bounds on the 
prices of barrier and exotic options. We have also shown how one can apply density estimation 
techniques for locally-elliptic random variables on the Wiener space and Malliavin calculus tools 
to prove that a lower bound of the same range holds the for the density of X as well. 

An central subject for future work is the way to generalise these results to a wider class of 
models, with particular focus on allowing for different powers of V in the equation for X and in 
the one for V. Our current impression is that this kind of program calls for a sharper version of 
the fundamental "tube" estimate for Ito processes which is at the basis of our analysis. 
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Appendix 



4.1 Some preliminary estimates 

Lemma 4.2. Assume (G) and let (X t , Vi;0 < t < T) be two processes satisfying f2.1]) - [£Jfy . 
Then for every < t < s < T and every p > 1 there exist a positive constant C p such that 



E 



sup (\X r - X t \ 2 P + \V r - V t \ 2p ) F t ] < C p (s - t)P exp(C p s 2 P) 



(4.20) 



Cp also depends on the parameters K,rj,a given in (G) and on Vq. 



Proof. Observing that both the functions v — > /3(t, v) and v — > o~(t, v)y/v have sub-linear growth 
under (G), (|4.20[) follows from the application of Burkholdcr's inequality and Gronwall's Lemma 
to the process (V t ; < t < T) satisfying (|2~T|) . then to (X t ; < t < T) satisfying 



4.2 Proof of Lemma 14.11 



Proof. Estimate (|4.4j) : Consider (x\,v\) G B[i 1 (y, \y\ + Vq)- On the set 

A^ s ' Vl (x xl ' y ,v Vl ' lyl+Vo ,R), with R u = R 2 ,ue [t,s], we have V u > \y\ -R 1 -R 2 > \y\ - 2y/\y\ 
and V u < \y\ + Ri + R 2 < \y\ + 2 ^/jiTT for all u G [t,s\. Therefore, ±\y\ < V u < 2\y\ for all 
u G [t,s], if \y\ > 16. Using estimates (fXTU]) . ([______]), and (fXl4) it is easy to show that 

conditions (|3.2|) . (|3.3[) and (|3.4[) are satisfied for the process (X*' X1 , V*' V1 ; t < u< s), the curves 

xi,y vi,\v\+Vo 



,V. 



\ the radius R. and the constant curves 



c\y\; 
: 7|y| 



A„ 



: L T \y\; 

pIMv\, 



(4.21) 



defined for u G [t, s], where c, Lt,7, A are the same as in Proposition 13.11 The derivatives a/,?/ 
and the radius R. being constant, all the involved curves belong to L(l,oo). The factor <f>\ tl in 

_____ 2 \ 2 

(|3.7p is still given by </>a. 7 = , hence Qa.7 = 1 ^2 hi with (7 = 8 12 e 2 . On the other hand, 
using the constant Ft defined in p.21j) . we have 

-«) 2 + K) 2 



F x , v ,n.(u)du = 
t Jt 



A 1, 



2« 



<4r T y_ ( 
< sr T ( 



R\ 



.{s-tfy 
R\ 



IT " 



1 



M 1 



!_)(«-*) 



( s - *)_, i?: 

Estimate Q3.5p then tells that ()4.4[) holds with the same constant ct defined in (|2.7[l . 

(|4~5|) : We fix i > 0, w G B eaV ^ /2 (y) and write V = V*'* for simplicity. We set B = {\V S - y\ < 

e VSy,t- 5 < s < t} and C = {|j____(<r(u,K) - *(* ~ <5> V._.))VKdW_ | < e V^}. Since 

e \/?y < £o\/So < 1, on the set B we have V. < 1 + y < 2y and \(3(s, V s )\ < K(2 + y) < 2Ky for 
all s G [t — 6, t] . Hence, on the set B 



\V s -v\ < 2KyS- 



t-s 



o(u,V u )^/V u N2ydWl 



and 2KyS < f 



> 



sup 

t-5<s<t 



Sy because S < So/y. Therefore 

V s - v\ < ^y/Syj > P( sup f a(u,V u )^/V u A2ydW^ < ^y/Sy). 



We time-change the stochastic integral into 6a 4) where (& s ; s > 0) is a standard Brownian motion 
and we denote A s — j t _ s o'(u,V u ) 2 (V u A 2y)du the quadratic variation (Dubins & Schwartz 
theorem, cf. Th. 3.4.6 in [T|)]). Since A s , s G [t — 6,t], is uniformly bounded by W 2 2yS, we 
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have sup t _ (5<s<t |6 J 4 s | < sup 0<s<25 :2 y5 |6 S |. Using the scaling propertyof the Brownian motion 
(b cs ; s > 0) ~ (^/cb s ; s > 0) we obtain 



P(B)>p( sup \b s \<^-^ 

K 0<s<2a 2 yS 4 



sup |6 S | < 

0<s<l 4v2cr 



3- 



The same arguments lead to 

P(SnC c )<p( / (a(u, V u ) - a(t - 5, V t -s))\/v~dWl 

V Jt-8 

>(|y (<r(u, K) - <r(t - S, V t - S )WVu A 2ydWl | > eov 7 ^) 
'(|6 Bt | >e ^) 



> e V^y; V u < 2y,t - S < u <t 



< 



where B t = f t _g((r(u, V u ) — a(t — 5, V t -s)) 2 {V u A 2y)du and (b s ; s > 0) is a standard Brownian 
motion. Using (R), we have B t < 2K 2 (5 2 + \V U - V t s\ 2 ) • 2y ■ 8 < AK 2 {8 2 + 4e 2 8y) ■ yS < 
2QK 2 4{8y) 2 1 hence 

¥(B n C c ) < 



sup |6 S | > e y^Sy 

s<20K 2 e 2 (Sy) 2 

(sup \b s \ > . 



< 20K 2 Sy E sup \b s \ 2 < 80K 2 S < -q 
L s <i J 2 

and we have used Doob's inequality and the value of So to get the two last inequalities. We 
conclude that 

P(B n C) = P{B) - P(C C nB)>q--q=-q. 

□ 



4.3 Conditional Malliavin calculus 

We briefly introduce the main elements of conditional Malliavin calculus. We consider a probabil- 
ity space (f2, F, P) with a filtration (Ft, t > 0), and a <i-dimensional Brownian motion (Wt, t > 0) 
with respect to this filtration. We fix some t > and S > 0. Conditional Malliavin calculus 
amounts to consider the standard Malliavin derivative operators, but focusing on the deriva- 
tives with respect to Wt,t £ [t, t + 5] on the one hand, and to replace expectations with con- 
ditional expectations with respect to Ft on the other hand. We recall the basic notation of 
Malliavin calculus (we refer to [35] for a more detailed presentation of this topic). D fe,p denotes 
the space of the random variables which are k times diffcrcntiable in Malliavin's sense in LP . Let 
0/c = {1, . . . , d} k be the set of multi- indexes of length k with components in {1, • ■ • , d}, and let 
jg>e fc _ {(x a ) aeQk : x a £ R}. For a measurable function V : [0,oo) fc — > W &k , we define 

\V\ 2 k := [ y \V a ( Sl ,...,s k )\ 2 d Sl ---ds k and 

H k := {V : [0,oo) fe R &k : \V\ 2 k < oo}. 

For F £ D fc,p and a £ <d k , the derivative of F of order k and index a is D k > a F. We denote 
D a u ...,s„ F = ( D au.-,s„ F )ceei. : il is known (cf. HI) that E[\D k F\ p k ] < oo for all p > 1, hence 

The conditional Malliavin calculus is based on the use of the following scalar product and 
norm: for every fixed t, S > and U, V : [t, t + S) -> R Bk we define 

(U, V) t ,6,k ■= f Y\ V a (s u . . . , s k )U a (s u . . . , s k )d Sl ■ ■ ■ ds k , 

\V\ls, k := (V, V) 2 tAk = [ Y, • • • • Sfc )! 2dsi ' ' ' dsfc - 

J[t,t+8)» ae&k 
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For F £ D fe,p , we define the following Sobolcv norms: 

El^M.^EE / \D i s '».. tSi F\ 2 ds 1 ---ds i , 



I FN 2 



i=o aee t J [t,t+sy 

l-^llt^.fcl-^]; 



i jpw\p 



Notice that H-FHt^j^p is not a constant (as happens in the standard Malliavin calulus) but an 
J-"t-measurable random variable. The standard norm corresponds to ||-F||o,oo,fc,p- Remark 

that by the definition of |||-F|||t,,5,fc,p, using the elementary inequality ^ < {a 2 + b)P/ 2 < 2P/ 2 6 p / 2 
for positive a, b and p > 2, we have 



i/<|||^|||^ p <2^ 2 /, 



(4.22) 



with / = E 



p/2 



"(eJU E Q£ e, I [t , t+S y \Dit.., Si F?dsi ■ ■ ■ 
We will make use of the two following inequalities: first, let F,G £ n p >iD fc,p . Then, for every 
P> 1, 

IMtAfc.P < W2 fc ||J ;, ||t,4 1 fc 1 2p||G||t, S)fc) 2p. (4.23) 

In addition, for every fc > 1 there exists a constant /x(fc) such that for every (j> £ C£, every p > 1 
and every F G D fc,p , one has 



\U(F)\\\ t , 5 , k , p </i(fc)|^U||| J F||| iAfe , 2fcp , 



(4.24) 



where |(/>|fc = X)i=o su PzeR l^^( x )l- Inequality (|4.24p is a consequence of the chain rule and of 
(|4.23p . The proof of (|4.23|) is based on some rather standard (but cumbersome) computations 
and can be found in [2J, Lemma 2.5 in the Appendix. 

Let us now consider diffusion processes. We consider some T > and < 5 < 1 AT and 
(Y t ; t £ [T — (5, T]) the unique strong solution to the equation 

Y t =Y T - S + B(s,Y s )ds + V / AjfaYJdWi, T-S<t<T, (4.25) 
Jt-s . =1 JT-8 

where y T _,5 G L 2 (ft, F T -s\ R") and B,Aj G Lip([T ~ S,T] x R n ;R n ) nC£°(R";lR™) for all 
7 = 1, . . . , d. We define 

e p (S) := e^'O^+WO*, 

and 

(fe+i) 2 



iv fc (A J B) = |^_ 1 gs| fe + |^| 

The following proposition gives the conditional version of the estimates in [BJ, Lemma 2.1 and 
Corollary 1. 

Proposition 4.4. For any k > 1 and any p > 1 there exists a positive constant d k , p depending 
on k,p but not on the bounds on B and A and their derivatives such that, for every T — 5 < t < T 
and for every I = 1, . . . , m, 



sup sup 

aee fc s u ...,s k £[T~8,T) k 



I r-)k,a -w-l \P 
\- u s 1 ,...,s k I t | 



l/p 



<d Kp N k {A 1 B)e p {8) d ^ 



\Y t l \\\T-s,8,k, P < 2kd k d Kp xVSx N k {A,B)e p {5) d ^. 



(4.26) 



(4.27) 



Proof. Inequality ()4.26() relies on the same proof as Lemma 2.1 in [BJ. Estimate (|4.27j) is a 
consequence of (|4.26|) : for any i = 1, . . . , k and I = 1, . . . , n, we have 

(/ \D^_ s jl\ 2 d Sl ---d Sl ) p/2 <5«%-V f irs_ y,' ",/>,...,/>,. 

X J\T-8.TY ' J\T-S.TY 
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hence, using (|4.22[) and Yli=i card(0;) = J2i=i d l — kd k , wc obtain 

\\\Y}\\\ P T- S ,5, k , P < 2 p/2 (kd k T . max max( / EO^"...^'^^^ ■ ■ • ds, 

f i=l,..., k a ^ e i^J[T-S,TY 

< 2 p ' 2 {kd k ) p x S p / 2 x d k , p N k (A,B) p e p (S) pdk ^, 
which proves (|4.27[) . □ 

4.4 Proof of Proposition 14.31 

Proof. We assume without further mention that we are on the set As^. Since T — 8 and 5 are 
fixed, we drop them from the notation and write \\F \\k, P instead of ||-F||T-5,5,fe,p and so on. 

We first show that estimate (|¥7T7|) holds for E[\R a \ p \Tt-8} 1/p ■ Using \r)(t, x)-r}(s, y)\\/ip(v) < 
Ky/y(js — t\ + \y — x\) and applying Burkholdcr's inequality we obtain 

i-T 



E[\R\*\Ft-s] < cM^ 1 f EMVtWFT-sW 
K Jt-s 

+ 5 p ' 2 - 1 f E[( V (t,X t ) - n(T - S,X T _ s )) p i,(V t ) p / 2 \^T-8}dt) 

JT-8 ' 

< c p ({5y) p + 5 p ' 2 - 1 K p y p ' 2 £ (5 P + \\X t - X T - S \\^ p )dt) 

< c p {8y) p ' 2 ({5y) p/2 + K P S P + K p C p 5 p / 2 e c ^ J 



<c p {8y) p e c ^\ 

where we have used 6 < \8 < ^ffry and Lemma [4.21 to estimate \\X t — Xt-s\% p = E[\X t — 
Xt-s\ p \^t~s\- Then we have 

||i?A||o. P <c p — ^e c " TP . (4.28) 
P±V 

We now estimate the Sobolev norms of X and V. Notice that by the definition of ijj we have 
\ip\o < U + § < 2y hence IVKOU — C ^V f° r au ^ — U f° r some constant c^. Similarly, Iv^lo < 



y + § < 2Vy and £ y/flv) = ^gL, hence 



< lV y lo < l^lo 



ish^l ™ 



and it can be seen easily that |v^?(')U — c< k ' ^ or an eventually different constant . With 
a slight abuse of notation, we write \a\fip\k (resp. I^lfc) f° r the | • |fc-norm of the function 
(t,v) — > a(t,v)y/tj}(v) (resp. (t, x, v) — > rj(t,x)^/ijj(v)). Then, using (|4.27[) with k = 2, for any 
t G [T - 5, T] we have 

|||Vt|||2,8p < 4 x d 2fip xVSx |(Tv^|?(|j8|a + WVi>hf 
xexp(d 2 , 8p 8**{\P\i + WV^P\i) ap ) 

< c^K 9 xVSx |v^l5lv^li x exp(d 2 , 8p S 4p K 16p \^P\l P ) 

< c {2) x y/Sy™ x exp( C ( 2) x (Sy) 4p ). 



Hence, using Sy < 1 wc get |||Vt|||2,8p < c p x Sy 11 , for a (eventually different) constant 

(2) 

c p '. An analogous estimate holds for |||X 4 |||2,8p: observing that the only difference is in the 
contribution of the drift term \rpp\2 < Kc$y, we have 

\\\X t \\\2,8 P < i 2) Vd X X |V^|? X CX P ( C ( 2 ) X (V^IlD 
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since \/~5y < 1, too. Now using (|4.24[) and denoting /i = p(2), we have ||??(t,^t) 2 ||2,2p < 

MMilll^tllksp and ||^(^)||2,2p < HV'blll^tllksp < A*4 2/|||Vt|||2,8p- Hence, using (|4.23|l . we 
get 

< 2!2Vi^ 2 cl 1) y x ll|X t ||| 2 . 



sup |||r7(f,X 4 )^(^ t )|||2, P < 2!2V^ 2 4 ; 2/ 

T-<5<t<T 



|V t ||| 2 ,8p ^CpV^V 1 , (4-29) 



where the constant c p also depends on K. We denote R = — |l+J, setting / = J^_ g r](t, X t ) 2 ip(Vt)dt 



and J = J T _ s {vt - ri T -s)Jip{V t ){pdWl + p±dW?). Then, using @22| 



IT-S\ 

T 2 

|/|||^<2^ 2 E^- 1 /(EE 

|||»?(t,X t )V(Ft)|||f )P di 



\P/2 




) dt 


F T -s 



T-S 



We now estimate the first Sobolev norm of J. For the ease of notation, we write rj t for r](t,X t )- 



For T - S < s < t < T, we have £ s lj ((?7 t - VT-s)\h{V t )) = ^'(Wl^t)), hence 



/ fat - nr- S )Jil>(Vt)(pdW} + P ±dW 2 ) = ( Vs - VT ^)M(V s ))p 3 

'T-S 



+ 



T-S 



D^( Vt ^(V t ))(pdWl + P± dW?) (4.30) 



with pj = plj=i + p±lj=2- Using bound ()4.26[) and proceeding as for (|4.29p we obtain 

i/p 



su Pi6[T-<5,T] 



E 



Ft -i 



< c p | r\ 1 1 1 i/0 1 i \J y 11 + 20 < c p \J j/ 32 , where c p depends 



also on K. Then, starting from (|4. 22[) and using Burkholder's inequalities 



+ ^E 



3=1,2 



E 



^W^ t ))| p 



dtds 



< VK vpl*-itfl* f (6+\\X a - X T _s\\o, P ) p ds + c p 6 p+1 y 32p/2 

JT-8 

<c p (5 2 y 32 ) p/2 e C * TP . 

The Sobolev norms of higher order are estimated in a similar way, giving the bound ||| J|||2,p < 
c p y/8 2 y 32 e CpTP . Finally, |||i?||| 2lP < |||/|||2,p + |||</|||2,p < c p \/S 2 y 32 e CpTP and this last estimate 
together with (|4~2"g|) yields |||-Ra|||2, p < c P j^V% 3 ~ I e CpTP , which is (j4~T7| . □ 
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